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I. INTRODUCTION 



In the past few years there has been considerable work on non-Hermitian Hamiltonians. 
Among this kind of Hamiltonians, much attention has been focused on the investigation of 
properties of so-called VT — symmetric Hamiltonians. Following the early studies of Bender 
et al. [1], the VT — symmetry formulation has been successfuly utilized by many authors 
[2-8]. The VT — symmetric but non-Hermitian Hamiltonians have real spectra whether the 
Hamiltonians are Hermitian or not. Non-Hermitian Hamiltonians with real or complex 
spectra have also been analyzed by using different methods [3-6,9]. Non-Hermitian but 
VT— symmetric models have applications in different fields, such as optics [10], nuclear 
physics [11], condensed matter [12], quantum field theory [13] and population biology [14]. 

Exact solution of Schrodinger equation for central potentials has generated much interest 
in recent years. So far, some of these potentials are the parabolic type potential [15], the 
Eckart potential [16,17], the Fermi-step potential [16,17], the Rosen- Morse potential [18], 
the Ginocchio barrier [19], the Scarf barriers [20], the Morse potential [21] and a potential 
which interpolates between Morse and Eckart barriers [22]. Many authors have studied 
on exponential type potentials [23,24,25,26] and quasi exactly solvable quadratic potentials 
[27,28,29]. In addition, Schrodinger, Dirac, Klein-Gordon, and Duffin-Kemmer-Petiau equa- 
tions for a Coulomb type potential are solved by using different method [30,31,32,33,34]. 
The exact solutions for these models have been obtained analytically. 

Further, using the quantization of the boundary condition of the states at the origin, 
Znojil [35] studied another generalized Hulthen and other exponential potentials in non- 
relativistic and relativistic regions. Domingues-Adame [36] and Chetouani et al. [37] also 
studied relativistic bound states of the standard Hulthen potential. On the other hand, Rao 
and Kagali [38] investigated the relativistic bound states of the exponential-type screened 
Coulomb potential by means of the one-dimensional Klein-Gordon equation. However, it 
is well-known that for the exponential-type screened Coulomb potential there is no explicit 
form of the energy expression of bound-states for Schrodinger [39] , KG [38] and also Dirac 
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[16] equations. In a recent work [31], §im§ek and Egrifes have presented the bound-state 
solutions of one-dimensional (ID) Klein-Gordon equation for VT — symmetric potentials 
with real and complex forms of the generalized Hulthen potential. In a latter study [32], 
Egrifes and Sever investigated the bound-state solutions of the ID Dirac equation for real 
and complex forms of generalized Hulthen potential for VT — symmetric potentials with 
complex generalized Hulthen potential. In recent works, we have solved the ID Schrodinger 
equation with the VT — symmetric modified Hulthen and Woods-Saxon (WS) potentials for 
£ 7^ bound-state spectra and their corresponding wave functions [40,41] using the Nikiforov- 
Uvarov (NU) method [42]. In the latter case, we investigated the VT — symmetric property 
and the reality of the spectrum for different real and complex versions of the modified WS 
potentials [41]. 

In the present study we investigate the bound-states solutions of the ID spinless Salpeter 
(SS) equation for the real and complex forms of the generalized Hulthen potential and 
correponding eigen functions using the same method. Thus, the objective of this work 
is to deal with the VT — symmetric property and the existence of bound-states (i.e., the 
reality of energy spectrum) when we solve the spinless Salpeter (SS) equation for some real 
and complex potentials with standard generalized forms. In view of the VT — symmetric 
formulation, we shall apply the NU method to solve the s-wave SS equation. In this regard, 
it is possible to present the theory of special functions by starting from a differential equation 
which is based on solving the SS equation by reducing it into a generalized equation of 
hypergeometric form. We also seek to present exact bound states for a family of exponential- 
type potentials, i.e., generalized Hulthen potential which is reducible to the standard Hulthen 
potential, Woods-Saxon potential and exponential-type screened Coulomb potential. This 
family of potentials have been applied with success to a number of different fields of physical 
systems. 

The organization of the present work is as follows. After a brief introductory discussion 
of the NU method in Section II, we obtain the bound-state energy spectra for real and 
complex cases of generalized Hulthen potentials and their corresponding eigenfunctions in 
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Section III. Finally, we end with some results and conclusions in Section IV. 



II. THE NIKIFOROV-UVAROV METHOD 

In this section we outline the basic formulations of the method. The Schrodinger equation 
and other Schrodinger-type equations can be solved by using the Nikiforov-Uvarov (NU) 
method which is based on the solutions of general second-order linear differential equation 
with special orthogonal functions [42]. It is well known that for any given ID radial potential, 
the Schrodinger equation can be reduced to a generalized equation of hypergeometric type 
with an appropriate transformation and it can be written in the following form 

where a(s) and a(s) are polynomials, at most of second- degree, and f(s) is of a first-degree 
polynomial. To find particular solution of Eq.(l) we apply the method of separation of 
variables using the transformation 

^n(s) = (f) n (s)y n (s), (2) 

which reduces equation (1) into a hypergeometric-type equation 

*Wni*) + T(s)y' n (s) + Xy n (s) = 0, (3) 

whose polynomial solutions y n {s) of the hypergeometric type function are given by Rodrigues 
relation 

Vn{s) = ^fy^T > (" = 0, 1, 2, ...) (4) 

where B n is a normalizing constant and p{s) is the weight function satisfying the condition 
[42] 

± w(8) = I|£) m(s) , (5) 
where w(s) = a(s)p(s). On the other hand, the function (f)(s) satisfies the condition 



d . , tt(s) , , . 

= -y<t>{ s ), (6) 

as a(s) 



where the linear polynomial 7r(s) is given by 



\ 



'^mV_ 5(s)+Ms) , (7) 



from which the root k is the essential point in the calculation of 7r(s) is determined. Further, 
the parameter A required for this method is defined as 

\ = k + n'{s). (8) 

Further, in order to find the value of k, the discriminant under the square root is being set 
equal to zero and the resulting second-order polynomial has to be solved for its roots 
Thus, a new eigenvalue equation for the SE becomes 

X n + nr'(s) + n{n ~ l) a"{s) =0, (n = 0, 1, 2, ...) (9) 

where 

t(s) =t(s) + 27t(s), (10) 
and it must have a negative derivative. 



III. RADIAL SPINLESS SALPETER EQUATION 

The relativistic wave Salpeter equation [43] is constructed by considering the kinetic 
energies of the constituents and the interaction potential. In addition, the spinless Salpeter 
(SS) for the case of two particles with unequal masses mi and rri2, interacting by a spherically 
symmetric potential V(r) in the center-of-momentum system of the two particles is given 
by 

^/-A + mf + V(r) - M 

_i=l,2 



x (r) = 0, A = V 2 



(11) 
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where the kinetic terms involving the operation ^ — A + mf are nonlocal operators and 
x(r) = Yt jm (9,(j))R nt i(r) denotes the Salpeter's wave function. For heavy interacting parti- 
cles, the kinetic energy operators in Eq. (11) can be approximated, (cf. e.g., Jaczko and 
Durand [44], Ikhdair and Sever [45,46]), as 

j: v '-A +m |= mi+m2 ----- .... (12) 

where a = mim2 stands for the reduced mass and ri = a I — mi?il2 „ 2 is an introduced 
useful mass parameter. This SS-type equation retains its relativistic kinematics and is 
suitable for describing the spin-averaged spectrum of two bound particles of masses m x and 
m 2 . The Hamiltonian containing the relativistic corrections up to order (v 2 /c 2 ) is called a 
generalized Breit-Fermi Hamiltonian (cf. e.g., Lucha et al. [47]). Hence, the SS equation 
can be further written in the form (h — c — 1) [45,46] 

{"2^ " + V(r) } Rnl(r) = EnlRnliT) ' (13) 
where E n i = M n \ — m\ — m 2 refers to the Salpeter binding energy with M n i is the 
semirelativistic-bound-state masses. 1 Moreover, it is worthwhile to point out here that, 
to obtain a Schrodinger-type equation, as could be seen fron Eq.(13), the perturbation term 
can be treated by using the reduced Schrodinger equation [49] 

p 4 = 4/i 2 [E nl - V(r)] 2 , (14) 

with p 4 = A 2 , and consequently one would reduce Eq. (13) into a Schrodinger-type form 
[45,46] 

[A /j, 2 



2/i 2rf> 



Ki + V \r) - 2E nl V{r) + V(r) R nl (r) = E nl R nl (r). (15) 



1 This approximation is correct to 0(v 2 /c 2 ). The A 2 term in (13) should be properly treated as a 
perturbation by using trial wavefunctions [48] . 



6 



In addition, the three-dimensional (3-D) space operator in the spherical polar coordinates 
takes the form 

d 2 2 9 L 2 

A = ^ + -^-^> ( 16 ) 

with L 2 = I (I + 1) . Hence, after employing the following transformation 

Rni\r) = — - — , (17) 

we obtain [45] 

A = ^_ [([+1) A 2 = _ 21 (j + 1) 92 + 4*(* + l) g + *(* + l)[Z 2 + Z-q (lg) 

Thus, using Eqs (17) and (18) and after a lengthy algebra but straightforward, we can finally 
write Eq.(15) in the following 3D space 



^ ^ | *(* + !) ft 2 , w(r) W nl (r) 2 



Mr) = 0, (19) 



where 



W nl {r)=V{r)-E nh m = rf / iS = ( mi m 2 /i)/(m im2 - 3/i 2 ). (20) 

It is worthwhile to point out that the Schrodinger-type equation (19), is found to be same 
as the formula given by Durand and Durand [50]. The perturbation term, W n i(r) 2 ; that is, 
(v 2 /c 2 ) term in Eqs. (19) and (20) is significant only when it is small (i.e., W n i(r)/m <C 1). 
This condition is verified by the confining potentials used to describe the present system 
except near the color— Coulomb singularity at the origin, and for r — > oo (i.e., the wavefunc- 
tion vanishes at and oo). However, it is always being satisfied on the average as stated by 
Ref.[50]. 

In the present work, we shall study the SS equation with a family of exponential potentials 
which is called generalized Hulthen potential [51], one of the important molecular potentials, 
in the ID -vector form, 

V q (x) = -V i ^*_ ax , V = Ze 2 a, < x < oo (21) 
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with a denotes the screening (range) parameter, Vo denotes the coupling constant and q is 
the deformation parameter which is used to determine the shape of potential. We have to 
note that, for some specific q values this potential reduces to the well-known types: such as 
for q — 0, to the exponential potential, for q — 1 to the standard Hulthen potential, and 
for q = —1 to the Woods-Saxon potential [41]. Further, near the origin, it reduces into the 
shifted linear potential in the limit of very short range (i.e., a — > 0) [31,32] 

V q (x) « -A- + r^-^ax + 0(«V), (22) 
q — 1 (q — iy 

with a constant shift, V /(q — 1) It also approximates to the screened Coulomb effective 
potential for small ax (i.e., ax — > 0) as [34] 

C"M « + (23) 

The complex form of the potential in (21) is said to be a VT -Symmetric potential when 
(cf. Refs.[7,31,32,41]) 

[VT ,V(x)\ = 0, (24) 

i.e., the VT — symmetry condition for the given potential V(x) satisfies 

[V(-x)Y = V(x). (25) 

To calculate the energy eigenvalues and the corresponding eigenfunctions, the Hermi- 
tian real- valued Hulthen potential form given by Eq.(21) is substituted into the ID 
VT — symmetrical Hermitian Schrodinger-type equation (19) for t = case (i.e., s-wave 
states): 



d 2 j) nq {x) 2fi 
dx 2 h 2 



El V e~ a * V 2 e- 2a * V E n e 
E n + + — — + + 



2m {l-qe~ ax ) 2m(l - qe~ ax ) 2 m(l-qe~ ax ) 



i>nq{x) = 0, 

(26) 



where \i = to/2 and to = 2m for two identical interacting particles. Now, introducing a 
convenient dimensionless transformation, s(x) = e~ ax , satisfying the arbitrary boundary 
conditions, 0<x<oo— >l<s<0, reduces Eq.(26) to the form: 
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d 2 ip nq (s) ldip nq (s) 2/j, 



+ 



+ 



ds 2 s ds h 2 a 2 s 2 



E„ 



El V n s 



+ 



V 2 s 



2„2 



+ 



V E n s 



2m (1 — qs) 2m(l — qs) 2 m(l — qs) 



^ng(s) = 0, 
(27) 



with the dimensionless definitions given by 



2/x Vb-g n 
fi 2 o; 2 m 



( e 3 > 0) 



(28) 



and finally one can arrive at the simple hypergeometric equation given by 



«.(»> + ^^.w + |s2(ei " fV ' y3 :r )+ ^ i+f3+29e2) ' e2| ^w = o- 



s(l - gs) 



[s (1 - 



(29) 



Hence, comparing the last equation with the generalized hypergeometric type, Eq.(l), we 
obtain the associated polynomials as 

f(s) = 1 - qs, a(s) = s(l - qs), a(s) = s 2 {e\ - q 2 e 2 - ge 3 - qe x ) + s (e x + e 3 + 2ge 2 ) - e 2 . 



(30) 



When these polynomials are substituted into Eq.(7), with <r'(s) = 1 — 2qs, we obtain 



ir(s) = -f ± \ ^s 2 (q 2 + 4(g2 e 2 + ~ + ~ _ e |) _ Aqk) + 4s (fc _ ~ _ e3 _ 2ge 2) + 4e 2. 



(31) 



Further, the discriminant of the upper expression under the square root has to be set equal 
to zero. Therefore, it becomes 



A = 



k - ei - e 3 - 2qe 2 - e 2 q 2 + 4(gV + qe x + ge 3 - e 2 2 ) - Aqk 



= 0. 



(32) 



Solving Eq.(32) for the constant k, we obtain the double roots as A; +i _ = e\ + e 3 ± 6e, where 
b = \Jq 2 — 4e|. Thus, substituting these values for each fc into Eq.(31), we obtain 
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as 1 (2ge - b) s - 2e; for k + = e 1 + e 3 + be, 
< s ) = —2 ± 2] 

{2qe + b) s — 2e; for k- = e\ + e 3 — be. 
Hence, making the following choice for the polynomial 7r(s) as 

7r ( a ) = -^-I[(296 + 6) a -2c], 

for fc_ = e\ + e 3 — 6e, giving the function: 



(33) 



(34) 



r(s) = -g(2 + 2e + 6/g)s + (1 + 2e), 



(35) 



which has a negative derivative of the form r\s) = — g(2 + 2e + 6/g). Thus, from Eqs.(8)-(9) 
and Eqs.(34)-(35), we find 



A = -|(l + 2e)(l + 6/5) + (ei + e 3 ), 



(36) 



and 



A n = (1 +n + 2e + b/q) nq. 



(37) 



Therefore, after setting \ n = A and solving for e, in K — 1 units, we find the Salpeter exact 
binding energy spectra as 



E na = \ — — m] < 
nq \2q 



1 ± 



1 - 



2ma 



2^,2 



\ q {f q ~™) D 



> , a 



2// ' 



(38) 



where 



D = (C 2 + 4e 2 ) /q = q + q(2n + l) 2 + 2(2n + 1)6, C = b + q(2n + 1) 



(39) 



where b = \jq 2 — for equal mass case. For convenience, when m\ = m 2 , the upper 
expression (38) can be further rearranged as 



E na — \ — — 2m | { 

nq \2q 



1 ± 



x (2mV f 



i - (-M + 1 M-V 

V 2mgV J ~ 4 V 2mqrVo / 



(g-2m)' 



> , < n < oo, 



(40) 
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where 



qa 



+ qa(2n + l) 2 + 2(2n + l)^q 2 a 2 -V 2 



= * 2 + vl 



K 



= ^fa 2 q 2 -V 2 + qa(2n + l), q 2 > ( 



(41) 



Let us now find the corresponding wavefunctions. Applying the NU method, the hyperge- 
ometric function y n (s) is the polynomial solution of hypergeometric-type equation (3) and 
described with the weight function [42]. By substituting n(s) and a(s) in Eq.(6) and then 
solving the first-order differential equation, we find 



(b+q) 

<f>(s) = s e (l - qs) 2 " . 



(42) 



It is easy to find the other part of the wave function from the definition of the weight function 



p{s)=s 2 \l-qsfl\ 



which then substituted into the Rodrigues relation resulting in 



(43) 



y nq (s) = D nq s e (l-qs) 



-b/q 



d n 
ds r 



-,n+2e 



l-qs) 



n+b/q 



(44) 



where D nq is a normalizing constant. In the limit q — > 1, the polynomial solutions of y n (s) 
are expressed in terms of Jacobi Polynomials, which is one of the classical orthogonal poly- 
nomials, with weight function (43) in the closed interval [0, 1] , giving y n ,i(s) ~ Pj> 2e ' b \l — 2s) 
[52]. The radial wave function ip nq (s) is obtained from the Jacobi polynomials in Eq.(44) 
and 4>(s) in Eq.(42) for the s-wave functions could be determined as 



^nq(s) = N nq s^(l - qs) — — s M (1 - qs) 



-,n+2e 



\n+b/q 



N nq s%l-qs)^P^\l-2qs), 

(45) 



with s = e ax and N nq is a new normalization constant determine by 
o o 



1 



/ \^ nq (s)\ 2 ds = N 2 q J s^l-qs)^ [P^\l-2qs)] 2 ds, 



(46) 
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We now make use of the fact that the Jacobi polynomials can be explicitly written in two 
different ways [52]: 

P<r\z) = 2- £(-ir p ( H + p )( n+ _ U ) (1 - zf-v (1 + zf , (47) 
p=0 \ P ) \ n PJ 

pM (z) - nn + P+l) "(n y(n + p + v + r + l) fz-lV 
n {) n \T(n + p + is+l)^ \r) T(r + p+l) \ 2 J' 1 1 

where (™) = r!( "_l r) , = r(r+ i)rff r+1) ■ Using Eqs.(47)-(48), we obtain the explicit expressions 
for P$*M*){\ -2qs) 

p{2e,b/ q){l _ 2gs) = (_i)«r(n + 2e + l)r(n + ^ + 1) 

n ( — \\Pqn-p 

X S P Kn- P y.T( P + l + inn + 2e- P+ lf" {1 " ^ (49) 

B ( 2gS) "r(n + 2 e + ^ + l) r ^ r!(n-r)!r(2e + r + l) (5 ° j 
Therefore, substituting Eqs.(49) and (50) into Eq.(46), it gives 
1 _ N 2 , ir+1 r (^+^ + l)r(n + 2 e + l) 2 f n (-W 



x 



r(n + 2e + | + 1) p\(n - p)\F(p + \ + l)T{n + 2e- P +l) 

l) r gT(n + 2e+* + r + l) 



S i? mp/o I ~ , n \ln q (p,r), (51) 

JZJ r!(n — r)!r(2e + r + 1) 1 



where 



I ng (p, r)=J s n+2 * +r ~ p (l - qsf+^ds. (52) 
o 

Using the following integral representation of the hypergeometric function [53] 

_ s) 7o-a -i (1 _ qs) -, ods =a Fi( ^ ; ) r(ao)r( 7o -ao) i 
^ r (7o) 

[Re(7o) > Re(a ) > 0, |arg(l - g)| < tt] (53) 
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which gives 

1 

2 F l {a ,(3 : a + l;q)/a = J s 010 ' 1 ^ - qs^ds. (54) 

o 

The hypergeometric function 2 Fi(a , j3 : 7 ; 1) reduces into 

7-1 / n n ,x r(7o)r( 7o -op- gp) 

2^l(«0,P0 • 70, 1) = 77 TFT 7TT; 

[Re( 7o - a - A)) > 0, Re( 7o ) > Re(A,) > 0], (55) 
for g = 1. Setting cto = n + 2e + r — p + 1, (3q — —p — ^ — 1, and 7 o = «o + 1, one gets 

T ( 2^1(00,^: 7o;g) (n + 26 + r-p+l)!(p+^ + l)! 

nglM a {n + 2 e + r-p + l){n + 2e + r + \ + 2)\- W 

A. Real potentials 

Firstly, we consider the real case in Eq.(21), i.e., all parameters (Vo, q, a) are real: 
(i) For any given a the spectrum consists of real eigenstate spectra E n (V , q, a) depending 
on q. The sign of V does not affect the bound states. It is clear that while V — > 0, 
E n = —2m 1 ± \l 1 — ( v ' l Zl >Ci ) which is for the ground state (i.e., n — 0) tend to the value 



E rs —3.73m and for the first excited state (i.e., n — 1) takes the value E 1 = —2m, where 
we have used A c = — = — = — which is the compton wavelength of the Salpeter particles. 

L mc ma c ° tr r 

(ii) There exist bound states (real solution) in case if the condition ^ < q 2 is achieved, 
otherwise there are no bound-states. 



(iii) There exist bound states in case if the condition (2m Vb) 

, 2 



xa 1 



< 



2mqVo 4 \ 2mqVo j 

— 2m) £, with £ = k 2 + V 2 is achieved, otherwise there are no bound-states. 
Moreover, this condition which gives the critical coupling value turns to be 

n ~ i " V ° ~ ^ ~ V ° 2 ) ~ \ (57) 

where 
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X 2 = 2q 2 



^ + ^ - 2 m 
Q \ 2( 1 i 



AmVn 



\ (--2m) + — 



(58) 



i.e., there are only finitely many eigenstates. In order that at least one level might exist, its 
necessary that the inequality 



qa+y/q*a*-V$<y/x i -VZ, 



(59) 



is fulfilled 

For a more specific case q = —1, Eq.(21) is reduced into the shifted Woods-Saxon (WS) 
potential 

V 



V(x) = -Vq + 



1 + e- 



(60) 



with energy spectrum given by 



E n (V ,a) = -(^ + 2m 



1 ± 



N 



1 - 



(2mV y 



i + i- 2 -) + ! M- 

' \ 2mV I ~ 4 \ 2mV 



(f + 2m 



> , < n < oo, 



(61) 



where 



a + a(2n + l) 2 - 2(2n + l)Ja 2 - V c 



(62) 



In this case, for any given a, all the eigenstates E n < 0. 

(iv) For q — 0, the potential (21) reduces into the exponential form 

V(x) = -V e~ ax , 



(63) 



the energy expression (40) does not give an explicit form, i.e., the NU method fails to give 
energy expression to this type of exponential potential. It is noted that for this potential 
there is no explicit form of energy expression of bound states for Schrodinger [16], KG [31,38] 
and also Dirac [36] equations. 
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For q — 0, the generalized equation of hypergeometric type which is given by Eq.(29) 
becomes 

Ms) + + [s24 + 3{€ \ +€3) -% M = 0, (64) 

s s 

with 

f(s) — 1, a(s) = s, a(s) = s 2 e\ + s(e 1 + e 3 ) - e 2 , (65) 
and the corresponding 7r(s) is determined as 

7r(s) = ± < 



(66) 



ie 2 s + e; for A; + = ei + e 3 + 2ie 2 e, 
ie 2 s — e; for fc_ = e 1 + e 3 — 2ie 2 e, 

where i = y^l and e 2 = Following a procedure similar to the previous case, when 
7r(s) = — ie 2 s + e is chosen for fc = e\ + e 3 — 2ie 2 e, then 

r(s) = -2ie 2 s + (1 + 2e), (67) 

and 

A = ei + e 3 - ie 2 - 2ie 2 e, X n = 2rae 2 , n (s) = s e e~* €2S , (68) 

could be obtained. Substituting a(s) and r(s), together with A into Eq.(3) gives 

sy^(s) + [1 + 2e - 2ie 2 s] y' n (s) - [ie 2 + 2ie 2 e - (ei + e 3 )] y n (s) = 0. (69) 

The last equation can also be reduced to standard Whittaker differential equation [54]. Thus, 
the solutions vanishing at infinity and it can be written in terms of confluent hypergeometric 
function as follows: 

Vn{s) = 1 F 1 ^ + e + * ( %* £3) ; 1 + 2e; 2ie 2 s) . (70) 
So, the acceptable solution for the upper component is found to be 

Ms) = <f>(s)y n (s) = A l F l Q + e + Z' 1 ^ ; 1 + 2e; 2ie 2 ^ s*e^ s . (71) 
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On the other hand, to find the an expression for the exact energy spectrum, when A n = A, 
in h — 1 units, we obtain 



„ . , 2n + 1 2m 2 
E n (a) = i — - — a — 2m — i- 



< n < oo. 



(72) 



2 ~ (2// • l)n' 

^From the last equation we conclude that if and only if ia is real then E n is a real. Hence, 
in order to apply NU method to this type of exponential potential e 2 should be complex 
(imaginary). This leads to result that either V or a must be imaginary. Therefore, for 
purely imaginary a, i.e., a — > iai, it reads: 

2m 



E n {a) = -m 



(2n+l) 
2 H ai + 



< n < oo. 



(73) 



2m ~" ' (2n+ l)a 7 

In the nonrelativistic limit, there is no available bound-state energy solution for the expo 
nential potentials [16]. 



B. Complex potentials 

Let us consider the case where at least one of the potential parameters be complex: 
(I) If a is a complex parameter (a — > ia), the potential (21) becomes 

V 



VJx) 



— [q — cos(o;a;) + % %m.{ax)\ = V*(—x) 



q 2 — 2qcos(ax) + 

which is a VT— symmetric but non-Hermitian. It has real spectrum given by 



E n (V , ia, q) = ( ^ - 2m\ ( 



1 ± 



1 + 



(2mV f 



x ~ V 2mgV / 4 V 2mqV J 



(f q -2m)' 



with 



q = q 2 a 2 + q 2 a 2 (2n + l) 2 + 2qa(2n + l)^/q 2 a 2 + V 2 , 



(74) 



(75) 



if and only if — (2mV o y 



1 + (asfc) + I (2^w)1 < ? (| - H 2 • The corresponding 



radial wave function ip nq (s) for the s-wave could be determined as 
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i> m {s) = N nq s ie (l - qs)^P^\l - 2qs), (76) 



where c = \J q 2 + and s = e~ iax . 

The norm of the wavefunction of such a non- Hermit ian quantum mechanical system is 
redefined as [2,55] 

oo 

Jr nq (-s)Ms)ds = u, u = ±l. (77) 



v — 1 stands for "PT— symmetric phase whereas v — — 1 stands for PT— antisymmetric 
phase. Therefore making the necessary parameter replacements in Eqs.(51)-(52), we can 
obtain the normalization constant for the complex VT— symmetric generalized Hulthen 
potential given by Eq.(74). 

For the sake of comparing the relativistic and non-relativistic binding energies, we need 
to solve the ID Schrodinger equation for the complex generalized Hulthen potential. We 
set the convenient transformation s(x) = e ~ tax ,0<x<oo^l<s<0, to obtain 

«.(.) + jf^-M + + y - ? ''-\ m = 0, (78) 

[s — qs z ) (s — qs z ) 

for which 

f(s) = 1 — qs, a(s) = s — qs 2 , a(s) = s 2 (qei — q 2 ^ 2 ) + s(2qe 2 — e\) — e 2 , 



e 2 = -^E n , e 1 = £ 1 V (ci>0). (79) 
ft cr ft cr 

Moreover, it could be obtained 

r(s) = -g(3 + 2?)s + (1 + 2e), (80) 

if 7r(s) = — ?(1 + e)s + e is chosen for k- = —qe — 1\. Finally, the bound-state spectra in the 
non-relativistic limit could be found as 



E n (V , q, ia) 



Ti 2 \2iiV /h 2 + qa 2 (n+ l) 2 



8fiq 2 a 2 



(n + 1) 



2 



> 0, 0<n<oo. (81) 
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On the other hand, for the sake of comparison, the non-relativistic limit for the real potential 
(21) can be found directly from the last equation as 

h 2 a 2 



E n (V ,q,a) = 



8/x 



(n + 1) 



(n + 1) 



, < n < oo. 



where = ^k. The radial wave function in the current case becomes 

' qa z n 

^ ng (s) = N nq s\l - qs)Pf' 1 \l - 2qs), 



with s = e iax and N nq is a new normalization constant determined by 

(n + l)\T(n + 2e + l) 2 | ^ (_i)p g "-J> 



1 = <(-I) n+1 " 



E 



r(n + 2e + 2) \^ p\(n-p)\(p + l)\T(n + 2e-p+l) 



(82) 



(83) 



x 



V (-l) r g r r(n + 2e + r + 2) | j +27+ _ +2 
r\(n-r)\T(2e + r + l) )J U gSj dS ' 



where 
i 



(84) 



J s n+2e+r-p^ _ gs )P+2 rfg =2 Fi ( n + 2? + r _ p + X> _ p _ 2 ; n + 21+ T - P + 2] q)B(n + 2t + T 


(85) 

(II) Let two parameters; namely, V and q be complex parameters (i.e., V — > iVo, g — > 
then the potential transforms to the form 



VJx) = V 



2cosh 2 (ax) — sinh(2aa;) — 1 


— i [cosh(o;x) - 


sinh(a:r)] 


1 + g 2 


2 cosh 2 (arc) - 


- sinh(2o;a;) — 1 





(86) 



which is a VT— symmetric but non-Hermitian if (2mVo) 
£ — 2rnj , it may possesss real spectra as 



\2mV q) ^ 4 \2mqV ) 



< 



E n (tV , a, iq) = ( ^ - 2m"j <i 



1 ± 



i-(-M + ±M-V 



where £ is defined by Eq.(41). On the other hand, the corresponding radial wave functions 
ip nq {.s) for the s-wave could be determined as 
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^ nq (s) = N nq s*(l - iq s)^rP^\l - 2tqs), (87) 

with s = e~ ax and d = \jq 2 - The integral I nq (p,r) = } s n+2e+r - p (l - iqs) p+ l +1 ds is 
given by 

Inqip, r) =2 Fi(n + 2e + r - p + 1, -p - ^ - 1 : n + 2e + r - p + 2; ig)S(n + 2e + r - p + 1, 1). 

(88) 

(III) When all the parameters Vb, ct and g are complex parameters (i.e., Vb — > iVb, ct — > ia, 
q —> iq), we obtain 

Vb 7T 



W = 



— [g — sin(o;a;) — « cos(arr)] = V^*(— — x). 



(89) 



q 2 — 2q sin(o;a;) + 

This potential is a pseudo- Hermit ian potential [27,56] having a tc/2 phase difference with 
respect to the potential in case (I), it is also a VT— symmetric, rj = P-pseudo-Hermitian 
(i.e., PTV q {x)(PT)- 1 = V q {x), with P = r] : x -> £ -x and T : i -> -i) but non-Hermitian 
having real spectrum given by 



E n (tV ,ia,iq) = - 2m\ { 



1 ± 



\ 



1 + 



(2mV 



1 + f g 1 + 1 (s- 



(f q -2m)- 



(90) 



where 



f = q 2 a 2 + gV(2n + l) 2 - 2ga(2n + l)y/q 2 a 2 + Vf 



(91) 



On the other hand, the corresponding radial wave functions ip nq (s) for the s-wave could be 
determined as 



^ m {s) = N nq s ie (l - iqs)^P^i\l - 2iqs), 



(92) 



with s = e tax and c is defined after Eq.(76). The integral I nq (p,r) = j s n+2te+r p (l — 

o 

iqs) p+7 i +1 ds is given by 

Q 

Inq(p, r) —2 Fi(n + 2ie + r — p + 1, — p 1 : n + 2ie + r — p + 2; iq)B(n + 2ie + r — p + 1, 1). 

(93) 
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IV. RESULTS AND CONCLUSIONS 



We have seen that the s-wave Salpeter equation for the generalized Hulthen potential 
can be solved exactly. The relativistic bound-state energy spectrum and the corresponding 
wave functions for the Hulthen potential have been obtained by the NU method. Some 
interesting results including the VT — symmetric and pseudo-Hermitian versions of the gen- 
eralized Hulthen potential have also been discussed for the real bound-states. In addition, 
we have discussed the relation between the non-relativistic and relativistic solutions and the 
possibility of existence of bound states for complex parameters. Finally, we have shown the 
possibility to obtain relativistic bound-states of complex quantum mechanical formulations. 
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